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Abstract 

Non-small cell lung cancers (NSCLCs) cause high nnortality worldwide, and the cancer progression can be activated by 
several genetic events causing receptor dysregulation, including mutation or amplification. MicroRNAs are a group of small 
non-coding RNA molecules that function in gene silencing and have emerged as the fine-tuning regulators during cancer 
progression. IVliR-lSBa is known as a key regulator in skeletal and cardiac myogenesis, and it acts as a tumor suppressor in 
various cancers. This study demonstrates that miR-133a expression negatively correlates with cell invasiveness in both 
transformed normal bronchial epithelial cells and lung cancer cell lines. The oncogenic receptors in lung cancer cells, 
including insulin-like growth factor 1 receptor (IGF-IR), TGF-beta receptor type-1 (TGFBRl), and epidermal growth factor 
receptor (EGFR), are direct targets of miR-133a. l\/liR-133a can inhibit cell invasiveness and cell growth through suppressing 
the expressions of IGF-1 R, TGFBRl and EGFR, which then influences the downstream signaling in lung cancer cell lines. The 
cell invasive ability is suppressed in IGF-1 R- and TGFBRl -repressed cells and this phenomenon is mediated through AKT 
signaling in highly invasive cell lines. In addition, by using the in vivo animal model, we find that ectopically-expressing miR- 
133a dramatically suppresses the metastatic ability of lung cancer cells. Accordingly, patients with NSCLCs who have higher 
expression levels of miR-133a have longer survival rates compared with those who have lower miR-133a expression levels. 
In summary, we identified the tumor suppressor role of miR-1 33a in lung cancer outcome prognosis, and we demonstrated 
that it targets several membrane receptors, which generally produce an activating signaling network during the progression 
of lung cancer. 
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Introduction 

Lung cancer mortality, especially in non-small cell lung cancers 
(NSCLCs), remains the leading therapeutic issue in cancer 
treatment worldwide [1,2]. Due to the diversity of histology and 
genetic changes in NSCLC, new drugs offer an improved overall 
survival only to a small proportion of the patients while the 
majority can only be treated with palliative chemotherapy [3]. The 
strategy of targeted therapy, which is based on the driver 
manipulation and druggable target detection that can quickly 
and precisely distinguish the molecular differences between each 
individual patient and offer personalized treatment, is the most 
promising advancement in the management of NSCLC patients 
[2]. 

The dysregulation of receptor signaling pathways, such as 
receptor tyrosine kinases (RTKs) and transforming growth factor- 
beta receptors (TGFBRs), has been identified as the major cause 



for lung cancer formation [4,5]. The genomic amplification, 
mutation or rearrangement of these receptors may lead to 
persistent activation, causing the activation of cell survival 
signaling and cellular transformation in addition to increasing 
the invasive ability as observed in various lung tumors. In NSCLC, 
up to .50% in adenocarcinomas of Asian patients harbor activating 
mutations of EGFR [6], and approximately 20% of squamous cell 
carcinomas have fibroblast growth factor receptor 1 (FGFRl) 
amplifications [7]. IGF-IR is highly expressed in epithelial 
differentiated NSCLC tumors, and therefore, blocking IGF-IR 
signaling inhibits the proliferation and survival of cancer cells [8] . 
TGF-fi signaling is mediated by two specific cellular serine/ 
threonine kinase receptors, namely TGFBRl and TGFBR2, and it 
can induce EMT in tumor cells in association with the acquisition 
of motility and invasive properties [9]. 
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MicroRNAs (miRNAs), a group of non-coding RNAs, can 
repress the expression of multiple target genes through endoge- 
nous RNA interference machinery [10,11]. MiRNAs are dysreg- 
ulated in a variety of human cancers, and specific signatures of 
aberrantly expressed miRNAs harbor diagnostic, prognostic and 
therap(;utic implications. Increasing evidence suggests that miR- 
NAs confer delicate biological processes and provides robustness 
via regulation of target networks [12]. In malignant lung cancer, 
oncogenic miR-135b promotes cell invasion and metastasis by 
modulating K\TS2, [3-TrCP, NDR2 and LZTSl in the Hippo 
pathway [13]. Conversely, the tumor suppressor miR-486 directly 
targets the components of insulin growth factor (IGF) signaling, 
including IGF-1, IGF-IR, and p85a, in lung cancer [14]. These 
studies indicate that a specific mlRNA may regulate multiple 
targets that belong to a signaling pathway or that share similar 
biological functions. Therefore, miRNAs may become therapeutic 
targets in lung cancer treatment. 

MiR-133a, which belongs to the miR-133 family, was first 
identified as a muscle-specific miRNA, and it plays an important 
role in myoblast proliferation and differentiation during embry- 
onic muscle development. Dysregulated miR-133a results in 
cardiovascular diseases and skeletal muscle dysfunctions, thereby 
indicating that miR-133a is required for the functional mainte- 
nance of muscle cells [15]. Recentiy, an increasing number of 
reports have shown that miR-133a acts as a tumor suppressor in 
several cancers, including prostate cancer, bladder cancer and 
head and neck squamous cell carcinoma [16-18]. In addition, the 
expression level of miR-133a is significantly reduced in lung 
squamous cell carcinoma [19], implying that miR-133a may also 
mediate lung cancer progression. However, the detailed mecha- 
nism of how miR-133a regulates the development of NSCLC has 
not yet been thoroughly validated. 

In this study, we explored the functions of miR-133a and 
identified multiple oncogenic receptors as its downstream targets 
in lung cancer. We provided evidence that miR-133a may 
downregulate several membrane receptors that activate and 
participate in cross-talk through downstream signaling using the 
same mediators in cancer cells. Our results also confirmed that 
miR-133a can suppress lung cancer invasion, metastasis and 
proliferation. Most importandy, the miR-133a expression level is 
significandy associated with NSCLC chnical outcome, implying 
that miR-133a has potential to be a therapeutic and prognostic 
target for lung cancer treatment. 

Materials and Methods 

Cell culture and antibodies 

A549 (a human adenocarcinoma alveolar basal epithehal cell 
Une) and HI 299 (a human non-small cell lung carcinoma cell line 
derived from the lymph nodes) cells were purchased from the 
American T}pe Culture Collection (ATCC, Manassas, Virginia). 
EKVX and H23 lung cancer cell lines were purchased from the 
Developmental Therapeutics Program of the National Cancer 
Institute (NCI, Bethesda, Maryland). The H441 cell line was a 
kind gift from Dr. Win-Ping Deng (Institute of Biomedical 
Materials and Engineering, Taipei Medical University) [20] . The 
BEAS-2B cell line was a kind gift from Dr. Reen Wu (Departments 
of Internal Medicine, University of California at Davis) [2 1] . The 
human lung adenocarcinoma cell lines CLl-1, CLl-1, CLl-5, and 
CL1-5-F4 were derived from in vitro transweU and in vivo metastasis 
selection as previously described [22]. The BES-2B, EKVX, H23, 
H441 and CL series cell lines were maintained in RPMI media 
supplemented with 10% fetal bovine serum (FBS). The H1299 and 
A549 cells were maintained in DMEM with 10% FBS. 



The primary antibodies used for immunoblot analysis were 
rabbit anti-EGFR (Santa Cruz Biotechnology, CA, USA), rabbit 
anti-IGF-lR[3 (Santa Cruz Biotechnology), rabbit anti-TGFBRl 
(Santa Cruz Biotechnolog) ), rabbit anti-FGFRl, rabbit anti-AKT, 
rabbit anti-phospho AKT (ser473) (Cell Signaling Technology, 
Inc., MA, USA) and mouse anti-P-actin (Santa Cruz Biotechnol- 
ogy) antibodies. 

Gene and mlRNA expression profiles 

The microRNA expression levels of lung adenocarcinoma were 
retrieved from the Cancer Genome Adas (TCGA) data portal. 
The level 3 data of miRNA-seq was downloaded. The miR-133a 
expression levels of 480 samples, which contained 434 tumors and 
46 normal lung tissues, were extracted. The log 10 magnitudes of 
reads per million (RPM) were applied to a t-test to estimate the 
statistical significance. To identify the miR-133a target genes, the 
GSE20028 and GSE26032 datasets, which contains expression 
profiles of miR-133a-transfected cancer cells vs. control cells, was 
downloaded and analyzed. The fold changes of EGFR, TGFBRl, 
FGFR and IGF-IR were estimated. 

Quantitative PGR analysis 

Total RNA was isolated using TRIZOL reagent (Invitrogen, 
Carlsbad, CA) according to the standard protocol. The mature 
miR-133a and endogenous control RNU48 were analyzed using 
TaqMan MicroRNA Assays (Applied Biosystems, Foster City, 
CA). Briefly, 30 ng of total RNA was reverse-transcribed via 
SuperScript-III Reverse Transcriptase (Invitrogen, Carlsbad, CA). 
The cDNA was amplified with a TaqMan 2 x Universal Master 
Mix (Applied Biosystems), and miRNA-specific real-time PGR was 
performed using an ABI 7500 real-time PGR system. 

Invasion assay 

TransweU chambers (8-|j,m pore size; BD Falcon, Franklin 
Lakes, NJ) were coated with the appropriate amount of 30 (Jg 
Matrigel (BD Biosciences, San Jose, CA), and 2.5x10* cells were 
suspended in 10% NuSerum-containing media (Gibco BRL, 
Grand Island, NY, USA), seeded in the chamber, and cultured for 
20 hours (CLl-5 and A549) or 24 hours (BEAS-2B and CLl-1). 
To prevent the activation of AKT, 1 xlO^ cells were pre-treated 
with LY294002 (50 |xM for the BEAS-2B cells) or AKT 1/2 
inhibitor (7.5 nM for the CLl-5 cells) (Sigma-Aldrich, St. Louis, 
MO) for 24 hours, and then cells were seeded in the chamber with 
LY294002- or AKT inhibitor-containing media. Cells that 
invaded the chamber from top to bottom were fixed with 
methanol and stained with a 50 |J.g/mL solution of propidium 
iodide (Sigma-Aldrich, St. Louis, MO). The propidium iodide- 
positive signal was quantified using th(- Analytical Imaging Station 
software package. Each sample was assayed in triplicate. 

Lentiviral vector transduction 

Pre-miR-133a-encoding sequences were subcloned into the 
pLKO-AS2.neo vector (obtained from the National RNAi Core 
Facility in Academia Sinica, Taipei, Taiwan), and lentiviral 
vectors were prepared in accordance with standard protocols. 
CLl-5 and A549 cells were infected by lentiviruses with the same 
multiplicity of infection (MOI = 3) in medium containing poly- 
brene (8 |j,g/mL). One day after infection, the cells were treated 
with G418 to derive a pool of neomycin-resistant clones. 

Luciferase reporter assay 

The 3'UTR regions of receptors from the cDNA of CLl-0 cell 
line were subcloned into pGL3 basic vector [3'UTR region of 
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receptors/number of miR-133a target sequence: EGFR (18-972/ 
2), FGFRl (147-625/2), IGF-IR (3852-4209/1), TGFBRl (1981- 
2287/1) and INSR (691-1099/1)]. One day before transfection, 
CLl-0 cells were seeded in 12-well plates at a concentration of 
2.5x10* per well. Next, 200 or 500 ng of the pLKO-AS2.neo 
vector or pLKO-AS2 miR-133a plasmid was co-transfected with 
50 ng of pGL3-target gene-3'UTR. The ReniUa luciferase 
plasmid (pRL-TK, Promega, Madison, WI) was co-transfected 
as a transfix'tion control. Cells were lyst'd 36 hours after 
transfection, and luciferase activity was measured using a Dual- 
Luciferase system (Promega, Madison, WI) according to the 
manufacturer's protocol. 

Human phospho-receptor tyrosine kinase array 

The phosphorylation levels of receptors were detected by the 
Human Phospho-Receptor Tyrosine Kinase Array Kit (R&D 
systems Inc, Minneapolis, MN) according to manufacturer's 
instructions. Briefly, capture and control antibodies were spotted 
in duplicates on nitrocellulose membranes and incubated over- 
night with 300 |a.g of protein lysate. The membrane was then 
washed extensively with the provided bufier and further incubated 
with a pan anti-phospho-tyrosine antibody conjugated to horse- 
radish peroxidase (HRP). After incubation, arrays were washed 
and visualized using enhanced chemUuminescence (ECL) (Amer- 
sham, Buckinghamshire, UK). 

Transfection and cell treatments 

Transient knockdown experiments were conducted with the 
human siRNA-SMARTpool of IGFIR, EGFR, TGFBRl and 
siControl (Thermo Scientific, Wilmington, DE) for oncogenic 
receptors; anti-miR inhibitor and anti-Cd (Invitrogen, Carlsbad, 
CA) for miR-133a. CLl-5 or A549 cells were transfected with 
indicated siRNAs or anti-miRNA inhibitor by RNAiFect trans- 
fection reagent (Qiagen, Valencia, CA) according to manufactur- 
er's protocols. The cells were used for invasion experiments 
48 hours post-transfection and cell numbers were counted 
72 hours post-transfection. 

BEAS-2B cells were seeded in 6-well plates at a concentration of 
1x10^ for treatment with 50 nM LY294002 (Sigma-Aldrich, St. 
Louis, MO). The cells were used for invasion and proliferation 
assays 48 and 72 hours post-treatment, respecti\'ely. 

Cells were seeded in 60-mm disiies at a concentration of 2x10' 
for treatment with rIGF-1 (R&D Systems Inc., MN, USA) or 
human TGF-|3-1 (PeproTech, Rocky Hill, NJ). After serum 
starvation for 8 hours, CLl-5 and A549 cells were treated with 
250 ng/ml rIGF-1 for 30 and 2.5 min, respectively. For TGF-P-1 
stimulation, A549 cells were treated with 0.2 ng/ml TGF-|3-1 for 
5 min. 

Experimental metastasis in vivo 

0.1 ml of PBS containing single-cell suspension of CL1-5/AS2- 
Neo vector (Ctl) or CLl-5/AS2-Neo-miR-133a-expressing cells 
(1 xlO**) was injected into the lateral tail veins of 6-week-old NOD 
SCID mice (LASCO). The death of tumor cell-injected mice were 
recorded, the remaining mice were sacrificed at 39 days, and their 
lungs were examined for metastasis (n = 9 per group). The lungs 
were removed and fixed in 10% formalin and the numbers of lung 
nodules were counted under a dissecting microscope. 

Ethics statement, clinical lung cancer samples and 
immunohistochemistry 

The Institutional Review Board of the Taichung Veterans 
General Hospital approved this investigation, and informed 



written consent was obtained from all patients involved in this 
study. Frozen lung cancer specimens from 112 patients who 
underwent surgical resection of NSCLC were analyzed for the 

expression of miR-133a. None of the patients had received 
adjuvant chemotherapy. MicroRNA expression profiling was 
performed using a TaqMan MicroRNA Assay Kit (Applied 
Biosystems, Foster City, CA) and the ABI PRISM 7900 Real- 
Time PGR System. MiR-133a expression was quantified in 
relation to the expression of small nuclear U6 RNA. 

Statistical analysis 

Data are presented as the mean ±SD. The differences between 
two groups were assessed using Student's ^-test, and the Kaplan- 
Meier method was used to estimate overall survival. Differences in 
survival between two groups were analyzed using the log-rank test. 
Multivariate Cox proportional hazard regression analysis with 
stepwise selection was used to evaluate independent prognostic 
factors associated with patient survival, and the expression of miR- 
133a, age, gender, tumor stage, and histology were used as 
covariates. All analyses were performed with SAS version 9.1 
software (SAS, SAS institute, Gary, NC). Two-tailed tests were 
used, and P values <0.05 were considered to indicate statistical 
significance. 

Results 

MiRlSSa is downregulated in malignant lung cancer cells 
and inhibits the cell invasive capacity 

In the TCGA lung adenocarcinoma database (https://tcga- 
data.nci.nih.gov/tcga/tcgaHome2jsp), we found significantiy re- 
duced expression levels of miR-133a in 434 tumor tissues 
compared with 46 normal lung tissues (Figure lA). To further 
evaluate the role of miR-133a during lung cancer progression, the 
expression level of miR-133a in transformed normal bronchial 
epitheUum cells (BEAS-2B cells) and five lung cancer cell lines with 
difierent invasive capacities was examined. We found that the 
expression levels of miR-133a were higher in BEAS-2B cells and in 
cells with low invasiveness (EKVX and H23) [23] than in those 
with high invasive ability (A549, H441 and HI 299) [24] 
(Figure IB). These data suggest that miR-133a is dysregulated in 
lung cancer and that the expression of miR-133a negatively 
correlates with cell invasiveness. 

To further examine whether miR- 1 33a is associated with cancer 
progression, we suppressed the expression of miR-133a by an anti- 
miR- 133a inhibitor in BEAS-2B cells and examined cell 
morphologv' and im'asiveness. Inhibition of miR- 133a expression 
increased the population of spindle-shaped cells and enhanced 
BEAS-2B invasive ability (Figure IC and ID), implying that 
reduced miR- 133a expression in transformed lung bronchial 
epithelial cells may enhance the invasive phenotype. 

We further examined miR- 133a expression using a panel of 
well-established lung cancer sub-lines, namely CLl-0, CLl-1, 
CLl-5 and CL1-5F4, with increasing invasive capacity [22]. We 
found that the expression of miR- 133a gradually decreased as the 
invasion ability increased (Figure IE). Next, we manipulated the 
expression of miR- 133a in high invasive cancer cell lines, namely 
CLl-5 and A549, and ectopic expression of miR-133a decreased 
their invasive capacities (Figure IF, right panel). Conversely, 
treatment with the miR- 133a inhibitor increased the number of 
invasive cells in CLl-1 cells (Figure IG, right panel). Thus, these 
results suggest that miR- 13 3a governs cell invasiveness in lung 
cancer cells. In addition, cell proliferation of CLl-5 and A549 
were decreased in miR-133a-expressing cells (Figure SIA), 
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Figure 1. The miR-133a tumor suppressor modulates cell invasive capacity in lung cancer. (A) The expression levels of mlR-133a In 434 
tumor tissues and 46 normal tissues from the TCGA lung adenocarcinoma database were examined. (B) The expression level of endogenous miR-1 33a 
in BEAS-2B, EKVX, H23, H441, A549 and HI 299 cell lines. IVIiR-133a expression level (C, left panel), as measured by real time RT-PCR, and cell 
morphology (C, right panel) in BEAS-2B cells were examined 72 hours after transfection of the anti-mlR-l 33a inhibitor. (D)The number of invasive cells 
was counted 24 hours after seeding cells that post-treated anti-miR-133a inhibitor (100 nM) for 48 hours in a transwell containing matrigel. (E) The 
endogenous level of miR-133a in CL1-0, CL1-1, CL1-5 and CL1-5F4 cells. (F) Measurement of the invasion in CL1-5 and A549 cells transiently infected 
with the AS2-Neo (Ctl) or AS2-Neo-miR-133a-expressing viruses (right panel). The expression level of miR-133a was normalized to that of RNU48, 
which was used as an internal control (left panel) (G) 48 hours after transient transfection of the anti-miR-133a inhibitor, the invasive capacity of CL1-1 
cells was determined (right panel). The expression level of miR-1 33a was normalized to that of RNU48, which was used as an internal control (left 
panel). 

doi:1 0.1 371/journal.pone.0096765.g001 
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suggesting that miR-133a may inhibit growth in these lung cancer 
cells. 

Multiple oncogenic receptors are directly regulated by 
miR-133a 

To investigate the mechanism of miR-133a-mediated suppres- 
sion of cancer cell invasion, we used TargetScan (http://www. 
targetscan.org/) to predict miR-133a downstream targets. Among 
them, a target network containing several receptor tyrosine kinases 
(RTKs), such as FGFRl, EGFR, IGF-IR, insulin receptor (INSR), 
and a serine/threonine kinase receptor, TGFBRl, are presented 
as possible candidates, all of which possess several putative miR- 
133a binding sites at their 3'- untranslated regions (3'-UTR) 
(Figure S2A). We further investigated the mRNA levels of these 
receptors in miR-133a-overexpressing cancer cell lines by com- 
bining the microarray database (http:/ /www.ncbi.nlm.nih.gov/ 
geo/query/acc.cgi?acc = GSE20028 and http://www.ncbi.nlm. 
nih.gov/geo/query/acc.cgi?acc = GSE26032) (Figure S2B). The 
mRNA levels were consistent with the in silico prediction results, 
suggesting that FGFRl, EGFR, IGF-IR and TGFBRl may be 
downregulated by miR-133a (Figure 2A). To evaluate whether 
miR-133a directly regulates the expression of these receptors, the 
3-'UTRs of these receptors containing the predicted miR-133a 
target sites were cloned into the luciferase reporter plasmid for a 
reporter assay. The luciferase activity of the constructs containing 
the 3'UTRs of EGFR, FGFRl, IGF-IR, TGFBRl were 
significandy reduced in the presence of miR-133a, although the 
fold change in FGFRl was not as dramatic as the others, and the 
luciferase activity of INSR was not affected by miR-133a 
(Figure 2B). We further evaluated whether TGFBRl may be 
dir('('th' regulated by miR-133a since this corr(;lati()ii has never 
been reported. The luciferase activity was decreased in the 
presence of miR-133a, but this repression was alleviated by the 
addition of an anti-miR-133a inhibitor (Figure S3A). These results 
indicate that miR-133a may regulate the expression of FGFRl, 
EGFR, IGF-IR, and TGFBRl. 

To test whether FGFRl, EGFR and IGF-lR-mediated 
signaling pathways are activated in CLl-5 and A549 cell lines, 
we detected their activating levels by phospho-RTK proteome 
assays and found that EGFR and IGF-IR, but not FGFRl, were 
activated in these cells (Figure S4A). To investigate whether miR- 
133a regulates the endogenous protein expressions of these 
receptors, we transduced CLl-.^ and A549 cells with miR-133a- 
expressing viruses and detected the protein levels of these 
receptors. As shown in Figure 2C, the expressions of EGFR and 
IGF-IR were decreased by the ectopic expression of miR-133a in 
CLl-5 and A549 cells (Figure 20). Downregulation of TGFBRl 
expression was obser\'ed in A549 cells but not in CLl-5 cells where 
the endogenous level of TGFBRl was low (Figure 2C). The 
FGFRl protein level was mildly downregulated, which corre- 
sponded to the luciferase results (Figure 2C). Moreover, the 
ec:topic expression of miR-133a decreased the IGF-IR activity in 
CLl-5 cells (Figure 2D). These results demonstrate that miR-133a 
may down-regulate die EGFR, IGF-IR and TGFBRl to different 
degrees in lung cancer cells. 

Knockdown of IGF-IR or TGFBRl inhibits lung cancer 
invasion and proliferation through AKT-mediated 
signaling 

Since we identified that miR-133a could suppress cell invasive 
and proliferative abilities, we further evaluated whether EGFR, 
IGF-IR and TGFBRl could modulate these functions in CLl-5 
and A549 cells by the receptors specific siRNA-sUencing approach 



(Figure 3A). Forty-eight hours after siRNA transfection, downreg- 
ulation of IGF-IR decreased the number of invading and 
proliferating cells in both CLl-5 and A549 cells (Figure 3A, 3B 
and 3C). In addition, TGFBRl -silenced A549 cells also reduced 
invasive ability and cell growth compared to the control cells 
(Figure 3B and 3C). Furthermore, knockdown of EGFR expres- 
sion induced a morphological change similar to the miR-133a- 
expressing phenotype (Figure S5A), and decreased cell prolifera- 
tion in CLl-5 cells (Figure 3C). However, EGFR knockdown did 
not influence cell invasiveness in neither cell lines (Figure 3B). 
These resuks indicate that IGF-IR and TGFBRl control cell 
invasiveness and cell proliferation while EGFR governs the latter 
solely, which raises the possibility that EGFR may activate a 
different downstream signaling pathway from IGF-lR/TGFBRl 
in CLl-5 and A549 cells. 

One of the main downstream signals of these receptors is the 
activation of AKT, which acts as a key regulator in cancer 
progression by promoting cell growth, anti-apoptotic cflFects, and 
cell invasion [25], and these regulations can be detected in CLl-5 
(Figure S6A) in A549 [26,27] ceU Hnes. To test whedier AKT- 
mediated signal is involved in the cell lines used in the present 
study, we first detected the phosphorylation status of AKT by 
knocking down the expression levels of each receptor (Figure 3D). 
Phosphor>'lated AKT was reduced when IGF-IR and TGFBRl 
were knocked down in CLl-5 and A549 cells (Figure 3D). 
Interestingly, inhibition of EGFR did not affect the AKT activity 
in these cells (Figure 3D). These results indicate that IGF-IR and 
TGFBRl may control cell invasi\'eness and cell proliferation 
through the activation of an AKT-mediat("d signaling. 

MiR-133a regulates AKT activity through suppressing 
downstream targets of IGF-IR and TGFBRl 

To evaluate whether miR-133a regulates AKT activity, we 
ectopically expressed miR-133a in CLl-5 and A549 cells. 
Overexpression of miR-133a decreased AKT phosphorylation in 
both CLl-5 and A549 cells (Figure 4A). We then suppressed miR- 
133a expression by using an anti-miR-133a inhibitor and found 
that inhibition of miR-133a enhanced AKT activation in BEAS- 
2B cells (Figure 4B). To further examine whether miR-133a 
suppresses invasion through AKT signaling, we treated BEAS-2B 
cells with the miR133a inhibitor and PI3K/AKT inhibitor, 
LY294002 (Figure S7A). Increased cell invasiveness by miR-133a 
inhibition was abrogated in the presence of LY294002 (Figure 4C). 
To validate whether miR-133a inhibits AKT signaling through 
repressing IGF-IR and TGFBRl, we next treated cells with their 
ligands. In the presence of IGF- 1 , the AKT phosphorylation level 
was reduced in miR-133a-overexpressing CL-1-5 and A549 cells 
(Figure 4D). In addition, TGF-|3-induced AKT phosphorylation 
was also reduced in the miR-133a-overexpressing A549 cells 
(Figure 4E). These results show that miR-I33a regulat(-s IGF-I- 
and TGF-P-mediated downstream signaling through downregu- 
lating their receptors. 

MiR-133a suppresses cancer metastasis in vivo 

To investigate whether the expression of miR-133a inhibits 
cancer metastasis in vivo, CLl-5/Cfl and CLl-5/miR-133a cells 

were inoculated directly into tlu' circulation of mice. Mice with 
CLl-5/Ctl cells developed more pulmonary' nodules than those 
injected with CLl-5/miR-133a cells (mean number of nodules: 
99.88 ±2 1.93 for line Cti and 39.11 ±16.76 for line miR-133a, p< 
0.01) (Figures 5A). Similarly, the mortality of mice with CLl-5/Cti 
cells was higher than that of mice with CLl-5/miR-133a cells, 
whereas all mice with miR-133a were stiU alive at 39 days after 
tumor injection (Figure 5B). Therefore, the pulmonary metastasis 
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and survival analysis of the CLl-5 murine model support the 
findings that miR-133a expression could suppress metastasis in vivo. 

MiR-l 33a acts as a prognostic indicator of the clinical 
outcome 

Because miR-133a suppresses multiple oncogenic receptors and 
inhibits lung cancer cell invasiveness, we explored the potential 
biological significance of miR-133a expression in lung cancer 
progression. We evaluated the correlation between the miR-133a 
expression profile and the overall survival in tumor specimens 
from 112 lung cancer patients. The characteristics of the 112 
patients are provided in Table SI, and miR-133a levels were 
measured by real-time RT-PCR. Kaplan-Meier analysis showed 
that low expression levels of miR-133a were significantly 
associated with poor overall survival (P = 0.0409; Figure 5C). 
Multivariable Cox regression analyses were used to evaluate the 
associations of various independent prognostic factors with 
patients' survival (Table 1). The results showed that the phenotype 
of high miR-133a expression was a protective factor when age. 



gender, and tumor histological type, were considered (hazard ratio: 
0.455, 95% CI: 0.233-0.888; p = 0.0209) (Table 1). This result 
suggested that patients who lia\(: tumors with higher miR-133a 
expression levels may have better clinical outcome. 

Discussion 

Metastasis is one of the most critical hallmarks of cancer, 
especially in lung cancer. There is an urgent need for the 
identification of metastatic-related factors and their underlying 
molecular mechanisms. MiRNAs play an important role in 
tumorigenic and metastatic progression, and they are involved in 
the pathogenesis and prognosis of NSCLC. Studies of prognostic 
markers in lung cancer have indicated that let-7-a and miR-221 
expression levels correlate with an improved survival [28] and that 
miR-34a and miR-155 expression levels correlate with lymphoid 
infiltration and a poor survival [29] . Here, we suggested that miR- 
133a may serve as a prognostic biomarker and that its expression 
level positively correlates with the overall survival in NSCLC 
patients. We have also found that miR-133a modulates cancer 
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Figure 3. IGF-1R and TGFBR1 mediate cell invasion through AKT signaling in lung cancer cell lines. (A) Measurement of the protein 
levels of receptors CL1-5 or A549 cell transfected with slEGFR, silGF-IR or siTGFBRI. (B) Measurement of the cell invasion ability of CL1-5 or A549 cells 
transfected with siEGFR, silGF-1 R or siTGFBRI for 48 hours. The number of invasive cells was counted 20 hours after the transfected cells were seeded 
and was presented relative to invasion of cells transfected with siCtl. (C) Measurement of the cell proliferation ability of CL1-5 or A549 cells 
transfected with siEGFR, silGF-1 R or siTGFBRI for 72 hours. (D) Representative Immunoblots show protein levels of pAKT (Ser473), AKT and p-actin in 
these cells. 
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invasiveness and the proliferation ability in lung cancer cells. We 
also identified that miR-133a can downregulate multiple targets 
and that these targets are all oncogenic membrane receptors, 
namely EGFR, IGF-IR, and TGFBRl (Figure 5D). The 
oncogenic roles of these receptors have been well identified, and 
they are all important therapeutic targets or molecular markers for 
cancer treatment [30,3 1] . 

MiR-133a is a key regulator of cardiomyogenesis and acts as a 
tumor suppressor in multiple cancer malignancies [32,33]. 
However, the molecular mechanism of miR-133a-involved lung 
cancer progression has not been clarified. Here, we found that 
miR-133a functions as a tumor suppressor and that its expression 
level is significantly downregulated in NSCLC cells, and this 
phenomenon was observed even in highly invasiveness cells. We 
found that miR-133a can suppress cancer cell invasion, metastasis 
and proliferation abilities, and we also found that inhibition of 
IGF-IR and TGFBRl reduces both cell invasion and cell 
proliferation. In addition, ectopic expression of miR-I33a in 
highly invasive cell lines decreases IGF-IR- and TGFBRl- 



mediated AKT activation, thereby implying that miR-133a may 
modulate cell invasion through downregulating IGF-IR- and 
TGFBRl -mediated AKT signaling. Although previous studies 
have indicated that miR-133a can inhibit cell invasion and 
proliferation through downregulating EGFR [34] , our in vitro lung 
cancer systems show that EGFR expression does not affect the 
invasion ability in CLl-5 and A549 cells. However, there is a miR- 
133a-induc;ed morphological change similar to the EGFR- 
knockdown phenotype (Figure S5A). In addition to a decrease in 
invasiveness, ectopic expression of miR-133a reduces the cell 
proliferation activity in CLl-5 and A549 cells (Figure SIA). In 
lung cancer, EGFR overexpression is associated with cancer 
proliferation [35] and downregulated EGFR causes cell growth 
inhibition in CLl-5 cells (Figure ?>C). These results suggested that 
miR-133a may modulate EGFR-mediated proliferation ability 
through regulating other signaling pathway. In summary, we 
propose that miR-133a can regulate oncogenic receptor-related 
cell invasiveness and proliferation in NSCLCs (Figure 5D). 
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IGF-IR and INSR belong to the IGF signaling pathway, which 
plays central roles in cell growth, differentiation, survival, cell 
transformation and cancer metastasis [8,36]. Unlike INSR 
expression, which contributes less to tumorigenesis, the malignan- 
cy of dy.sregulated IGF-IR in cancer progression has been 
investigated, and targeting IGF-IR in cancer treatment is 
currently under development. However, the results of recent 
clinical trials with IGF-IR inhibitors have encountered several 
hindrances, which need to be resolved. For instance, due to the 
homology between IGF-IR and INSR, the inhibition of IGF-IR 
signaling disrupts INSR signahng and results in unwanted side 
effects, such as hyperglycemia and hyperinsulinemia. In this study, 
we noticed that both INSR and IGF-IR contain the miR-133a 
target sequence in the 3'UTR. However, miR-133a specifically 
targets IGF-IR and reduces its receptor activity while INSR is not 
affected. Therefore, miR-133a mimics may be promising small 
compounds in regulating IGF-IR signaling in tumor cells. 

In this study, we identified TGFBRl as a novel downstream 
target of miR-133a. However, TGF-(3 signaling has dual functions 
in cancer progression [9]. TGF-P has a potent growth inhibitory 
effect on epithelial and lymphoid tissues at the early stage of tumor 
development, but it also becomes a pro-oncogenic factor that 
stimulates tumor cell growth and invasiveness at the later stage [9] . 



We observed that the A549/TGFBR1 -silenced cells present an 
epithelial-like morphology with reduced cell invasive and prolif- 
erative ability, suggesting that TGFBRl has the oncogenic 
property in highly invasive cells. Both miR-133a overexpression 
and TGFBRl silencing leads to AKT activation, and miR-133a- 
overexpressing cells have reduced TGF-p/TGFBRl -mediated 
AKT signahng, thereby indicating that miR-133a may have the 
potential role to regulate TGF-P-mediated EMT. 

Previous studies have indicated that miR-133b-regulated 
FGFRl -mediated tumor growth occurs in gastric cancer [37]. 
Although our results showed that functional FGFRl signahng is 
not significandy changed in our lung adenocarcinoma cell system, 
FGFRl suppression by miR-133a is detected in the reporter assay 
(Figure 2B), and miR-133a-suppressed FGFRl mRNA levels have 
been detected in four squamous cell carcinomas (Figure S2B). In 
NSCLC patients, FGFRl amplifications regulates cell prolifera- 
tion through activation of the MAPK and PI3K/AKT pathways 
in about 20% of squamous cell carcinomas but they are rarely 
detected in adenocarcinomas [4]. Therefore, although FGFRl 
mediated regulation has not been detected in CLl-5 and A549 
cells, it may be an important factor in mLR-133a-mediated tumor 
suppression in squamous cell carcinoma tumorigenesis. 
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Figure 5. iVliR-133a expression decreases cancer metastasis in vivo and associates with longer survival rates in NSCLC patients. (A) 

Effects of miR-1 33a in metastasis in vivo. Left panel: histological examinations of the lung tumors by hematoxylin-eosin (H&E) staining. Representative 
lungs of mice intravenously (i.v.) injected with CL1-5/Ctl, CL1-5/miR-133a. White arrows indicate the lung tumors. Scale bars, 200 |im. Right panel: 
Numbers of metastatic tumor nodules were calculated 39 weeks after tail-vein injection (n = 9 per group). (B) The number of death/total (top panel) 
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(Mantel-Cox) test. P = 0.01 64 (C) Kaplan-Meier plots of overall survival in 112 NSCLC patients in high- and low-risk groups based on miR-1 33a 
expression levels (P = 0.0409). (D) Model showing that miR-1 33a modulates the receptor-mediated cancer malignancy in NSCLC. 
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In summary, the microRNA, miR-1 33a, inhibits the expression 
levels of multiple oncogenic receptors, i.e., EGFR, IGFIR, and 
TGFBRl, and mediates cell growth and invasion in lung 
adenocarcinoma. Moreover, miR-1 3 3a expression levels predict 



the clinical outcome in NSCLC, thereby serving as a suitable 
prognostic factor. Therefore, miR-1 33a is a promising theranostic 
agent, and the development of miR-1 33a mimics should possess a 
great potential for cancer treatment. 



Table 1. 


Multivariable Cox regression analysis of nniR-133a levels and overal 


survival in 112 NSCLC patients. 




Variable 


Hazard ratio (95% C.I.) 


p-value 


miR-1 33a 


0.455 (0.233-0.888) 


0.0209 


Gender 


0.711 (0.307-1.647) 


0.426 


Age 


1.057 (1.018-1.097) 


0.0038 


Stage 


2.5 (1.678-3.726) 


<0.0001 



doi:l 0.1 371 /Journal.pone.0096765.t001 
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Supporting Information 

Figure SI Cell proliferation of miR-133a-overexpress- 
ing CLl-5 and A549 cells. (A) Measurement of cell prolifer- 
ation of miR-133a-overexpressing CLl-5 and A549 cells. 72 hours 
after transient infection with AS2-Neo (Ctl) or AS2-Neo-miR- 
133a-expressing viruses, the cells (1x10"^) were seeded in 6-weU 
plates, and the cell numbers were counted on the indicated days. 
Each conditions were performed in triplicate. 
(TIF) 

Figure S2 FGFR, EGFR, IGF-IR and TGFBRl are direct 
targets of miR-133a. (A) Putative miR-133a binding sites of the 
five receptors, including FGFR, EGFR, IGF-IR, INSR and 
TGFBRl, were identified by computational altjorithms from 
Targetscan (a-e). (B) Fold changes of FGFR, EGFR, IGF-IR and 
TGFBRl were estimated from expression profiles (Data set: 
GSE20028 and GSE26032) of miR-133a-transfected cancer cells 
vs. control cells. SAS and HSC3 were the oral squamous cell lines; 
HI 57 and PC 10 were the lung squamous cell carcinoma cell hnes. 
(TIF) 

Figure S3 TGFBRl 3'UTR is suppressed by miR-133a. 

(A) Co-transfection of CLl-5 cells with AS2-Neo vector (Ctl) or 
AS2-Neo-miR-I33a-expressing plasmid with firefly luciferase 
fused with 3'UTR of TGFBRl for 6 hours and then incubated 
with anti-miR-Cd or an anti-miR-133a inhibitor (100 nM) in 
complete medium for 36 hours. Luciferase activity was measured, 
and the relative ratio of the activity in the miR-133a groups to that 
in the control vector group is presented. 
(TIF) 

Figure S4 Phospho-receptor detection in CLl-5 and 
A549 cell lines. (A) CLl-5 and A549 cell lysates were incubated 
with nitrocellulose membranes that were conjugated with 
phospho-receptor antibodies in duplicate. The phosphorylation 
levels of EGFR, FGFR, INSR and IGF-IR were determined. 
(TIF) 

Figure S5 Cell morphology of miR-133a-overexpressing 
or oncogenic receptor-silenced CLl-5 or A549 cells. The 

cell morphology of CLl-5 (A) or A549 (B) was determined 
72 hours after transient infection with AS2-Neo (Cfl) or AS2-Neo- 
miR-133a-expressing viruses (upper panel) or 48 hours after 
transient transfection with siCtl, siEGFR, silGF-lR or siTGFBRl 
treatment (middle and bottom panel, respectively). 
(TIF) 
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were pre-treated with the AKT inhibitor (7.5 |a,M) for 24 hours, 

and then seeded into chambers with AKT-inhibitor containing 
medium. The invasive cells were determined 20 hours post- 
incubation. (B) The cell numbers of CLl-5 were counted after 
AKT inhibitor treatment for 3 days. 
(TIF) 

Figure S7 Up-regulation of phospho-AKT mediated by 
an anti-miR-133a inhibitor can be reduced by the PI3K/ 
AKT inhibitor in BEAS-2B cells. (A) Representative immu- 
noblots showing the protein levels of pAKT (Ser473), AKT and P- 
actin in BEAS-2B cells after treatment with an anti-miR-I33a 
inhibitor (100 nM) with eidier DMSO or LY294002 (50 \lM) for 
48 hours. 
(TIF) 

Table SI MiR-133a expression in relation to clinical 
parameters and pathological characteristics. The clinical 
characteristics of the 112 patients with NSCLC are summarized. 

(DOCX) 

Materials and Methods SI Luciferase reporter assay 
with anti-miR-133a treatment. One day before transfection, 
CLl-5 cells were seeded in 12-well plates at a concentration of 
6x10* per well. Next, 200 ng of the pLKO-AS2 neo vector or 
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measured using a Dual-Luciferase system (Promega, Madison, 
WI) according to the manufacturer's protocol. 
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